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Abstract — In this article new bounds on weighted p-norms of am- 
biguity functions and Wigner functions are derived. Such norms 
occur frequently in several areas of physics and engineering. In 
pulse optimization for Weyl-Heisenberg signaling in wide-sense 
stationary uncorrelated scattering channels for example it is a key 
step to find the optimal waveforms for a given scattering statistics 
which is a problem also well known in radar and sonar waveform 
optimizations. The same situation arises in quantum information 
processing and optical communication when optimizing pure 
quantum states for communicating in bosonic quantum channels, 
i.e. find optimal channel input states maximizing the pure state 
channel fidelity. Due to the non-convex nature of this problem 
the optimum and the maximizers itself are in general difficult 
find, numerically and analytically. Therefore upper bounds on 
the achievable performance are important which will be provided 
by this contribution. Based on a result due to E. Lieb [1], the 
main theorem states a new upper bound which is independent 
of the waveforms and becomes tight only for Gaussian weights 
and waveforms. A discussion of this particular important case, 
which tighten recent results on Gaussian quantum fidelity and 
coherent states, will be given. Another bound is presented for 
the case where scattering is determined only by some arbitrary 
region in phase space. 

I. Introduction 

Time-frequency representations are an important tool in signal 
analysis, physics and many other scientific areas. Among them 
are the Woodward cross ambiguity function A 97 (t, v), which 
can be defined as ( T denotes complex conjugate) 



A 97 (r, v) 



and the Wigner distribution W 97 (r, v) 

t 



W S7 (t,^) = / g( T + -)j( T - -)e 



'dt 



(1) 



(2) 



r (2t, 2v) where 



T ,x 2 

where the functions g,j : K — > C assumed to be in £2 
Both are related by W gi (r,i/) = 2 A, 
7~(£) = 7(— t). Hence all results which will presented later on 
apply on Wigner functions as well. Due to non-commutativity 
of the shifts in t and v (in phase space) there exists many 

'Which can be relaxed to other spaces by the Holder inequality 



definitions of these functions which differ only by phase 
factors. In considering norms only, the ambiguities due to these 
phase factors are not important. 

To be consistent with the previous work in [2], [3] in this 
article the alternative definition 



A g7 (x) ^ (g,S x j) = I g(t){S xl ){t)dt 



(3) 



is used, where S x is the time-frequency shift operator given 

as 

(4) 



and x = (x±,X2) S K 2 . 

Note that A gl (x) = e' l ' KXlX2 A- g j(— x\, — X2). These opera- 
tors establish up to phase factors an unitary representation 
of the Weyl-Heisenberg group on £2^) — the so called 
Schrodinger representation (see for example [4]). They equal 
(again up to phase factors) the so called Weyl operators 
(Glauber displacement operators), i.e. perform phase space 
displacements in one dimension. 

It is an important and in general unsolved (non-convex) 
problem in many fields of physics and engineering to find 
normalized function g and 7 such that the following integral 



\{g,S xl )\ 2 C{x)dx 



\A gi {x)\ 2 C{x)dx (5) 



is maximized where C(x) could be some probability distri- 
bution and dx is the Lebesgue measure on R 2 . For example 
in radar and sonar application (5) is typically related to the 
correlation response with some filter g of a transmitted pulse 7 
after passing through a non-stationary scattering environment 
characterized by some C( ). This formulation is obtained for 
so called Weyl-Heisenberg signaling in wide-sense stationary 
uncorrelated scattering (WSSUS) channels [5], [2], [6], [7] 
where C(-) is called the scattering function. 
If considering 7 as a probability wave function in quantum 
mechanics (5) can be considered also as its overlap with 
some wave function g after several phase space interactions. 



In quantum information processing (5) is typically written as 
pure state fidelity 

(5) = Tr{IT g J S x U y S* x C(x)dx} = Tr{n 9 A(n 7 )} (6) 

where IT/ = f |/||^ 2 (/, •)/ is the rank-one projector onto / 
and Tr(-) denotes the trace functional. The middle term in (6) 
is the Kraus representation [8] of a bosonic quantum channel 
^.(•) [9], [12] which maps the input state LT 7 (rank-one density 
operator) to the output A(II 7 ). Minimizing the probability of 
error P e = 1 - Tr{LI g ,4(II 7 )} (see for example [10]) for 
rank-one measurements is then the maximization of the pure 
state fidelity, i.e. the following optimization problem: 

maximize Tr{n ff A(II 7 )} (7) 

9,7 

For each 7 the operator A(II 7 ) is a positive semi-definite trace 
class (thus compact) operator, hence (7) likewise reads 

maximize A max (A(X)) (8) 

TrX=l,X>0 

where the rank-relaxation follows from convexity of the 
maximal eigenvalue A max (-) and linear-convexity of A(-) (see 
for example [2], [3]). 

In a slightly more general context this paper considers 
|||A g7 | r C||i which directly gives the wei ghted v— norms of 
ambiguity functions in the form of 

||A S7 || r , c = (J \A gi (x)\ r C(x)dx) 1/r = \\\A gi \ r C\\l /r 

(9) 

where C : R 2 — > R + is now some arbitrary weight function. 
For r = 2 the results match then the examples given so far. 
Note furthermore that this topic is also connected to Renyi 
entropies H(r) of time-frequency representations 

H(r) = -^log||A S7 C^ (10) 
1 — r 

i.e. a measure of time-frequency information content [11]. 
II. Main Results 

The results are organized in a main theorem presenting the 
general upper bound to |||A g7 | r C||i. Then, two special cases 
are investigated in more detail. The first is dedicated to the 
overall equality case in the main theorem and important 
for Gaussian bosonic quantum channels. The second case 
discusses another application relevant situation motivated by 
WSSUS pulse shaping in wireless communications. But before 
starting, the following definitions are needed. 

Definition 1 Let < p < 00. For functions f : R — > C and 
F : R 2 -► C 

H/|]/= (/ \f(t)\ p dt) " \F\v={j \F{x)?dx^ P 



are then the common notion of p-norms, where dt and dx are 
the Lebesgue measure on R and R 2 . Furthermore for p = 00 
is 

ll/IU ^'ess sup \f(t)\ \\F\U ^'ess sup \F(x)\ 

If || is finite f is said to be in C P (R) (similarly if \\F\\ p is 
finite, F is said to be in C P (M. 2 )). 

For discussion of the equality case for the presented bound 
the formulation "to be Gaussian" for functions / : R — > C 
and F : R 2 -> C is needed. 

Definition 2 Functions f(t) and F(x) are said to be "Gaus- 
sian" if for a,b,c,C £ C, A e <C 2x2 and B E C 2 

f(t) = e - at2+bt+c F(x) = e -<*,Ac>+(B,x)+o (n) 
and Re{a} > and A* A > 0. 

Two Gaussians f(t) and g(t) are called matched if they have 
the same parameter a. 

The main ingredient for the presented analysis is the following 
theorem due to E. Lieb [1] on (unweighted) norms of ambi- 
guity functions. 

Theorem 3 (E. Lieb) Let A g7 (x) = (g,S x j) be the cross 
ambiguity function between functions g G C a (R) and 7 G 
£ b (R) where l=i + i. If2<p<oo with q = ^ < a < p 
and q < b < p, then holds 

||A iry |||<H( P> o,6)||ff||S||7llS (12) 

where H(p,a,b) = c p q (c a / q c b/q c p / q ) p/q , c p = 
p 1 /( 2 p)q- 1 /( 2 i)_ Equality is achieved with g and 7 being 
Gaussian if and only if both a and b > p/(p — 1). In 
particular for a = b = 2 

l|A 97 ||£<^|| 5 07||i (B) 



Actually Lieb proved also the reversed inequality for 1 < 
p < 2. Furthermore, for the case p = 2 it is well known that 
equality holds in (13) for all g and 7. Then the optimal slope 
(related to entropy) 

-J \A gl {x)\ p \n\A g7 (x)\ p dx (14) 

at p = 2 is achieved by matched Gaussians [1]. For simplifi- 
cations it is assumed from now that \\gW2 = W1W2 = 1- With 
the previous preparations the main theorem in this article is 
now: 



Theorem 4 Let A S7 (x) = {g,S x "f) be the cross ambiguity 
function between functions 5,7 with ||g|| 2 = H7II2 = 1 and 
p,r6i Furthermore let C(-) £ C q (R 2 ) with q = TTien 



r C||i< 



(15) 



holds for each p > max{l, -}. 

Proof: In the first step Holder's inequality gives 

lllA^rCllx^lHA^nipllCII, (16) 

for conjugated indices p and q, thus with 1 = ~ + |. Equality 
is achieved for 1 < p < 00 if and only if there exists A £ K 
such that 



(17) 



for almost every .t. Similar conclusions for p = 1 and p = 00 
are not considered in this paper. Lieb's inequality in the form 
of (13) for ||A fl7 ||^ gives for rhs of (16) 



1 Airmen, 



|A ff7 ||; p ||C|| g 



= (||A fl7 p'||C|| g <(l)'||C||, 

(18) 

The latter holds for every rp > 2, thus the case rp = 2 is 
now included as already mentioned before. Equality in (18) 
is achieved if g and 7 are matched Gaussians. Furthermore if 
strictly rp > 2, equality in (18) is only achieved if g and 7 
are matched Gaussians. Replacing q = gives the desired 
result. ■ 

Note that apart from the normalization constraint the bound in 
Thm.4 does not depend anymore on g and 7. Hence for any 
given C( ) the optimal bound can be found by 



mm 

K9p>max{l,f} 



(19) 



In the minimization p > 1 has to be forced for Holder's 
inequality and p > - for Lieb's inequality. Two special cases 
are investigated now in more detail which are relevant for 
application. 

First the overall equality case in Thm.4 is considered. 

Corollary 5 Let C(x) = ae -°^(A+A) w ith R 3 a > 0. 
Then for each p > max{l, ^} holds 



2a \ p ( p — 1 



(20) 



The best bound is given as 

la 



a > 



2-r 



"Clli < { 2a+r ~ 2 

ai(l -r/2) 1 - r / 2 else 



(21) 



For a > — 2^ equality is achieved at p = — + 1 if and only 
if g and 7 are matched Gaussian, i.e. then 

2a 



2a + r 



(22) 



holds. 



Proof: The moments of L\ -normalized two-dimensional 
Gaussians are given as 



According to Thm.4 the upper bound 



del 



ICII 



2a\ p ( p — 1 



(23) 



(24) 



holds for each p > max{l,2/r}. The optimal (minimal) 
bound is attained as some point p mul which can be obtained 

as 

(25) 



min , f(p)= /(Pmin) 
max{ 1 , - } 



The first derivative /' of / at point p is 

ft, \ f(p) , / Kg" 1) n 

/ (P) = — 2~M 7T ) 

p z la 

Thus /'(p m i n ) = gives only one stationary point p„ 



1) 



1 <^ p 



2a 

mm - + 1 > 1 



(26) 



(27) 



2a r 
Due to f{p)/p 2 > and strict monotonicity of ln(-) follows 
easily that /'(p min + e) > > /'(p min - e) for all e > 0, 
hence / attains a minimum atp m i n . The constraint p m i n > 1 is 
strictly fulfilled for every allowed a and r, hence the solution 
is feasible (p m j n > -) if a > Then the optimal (minimal) 
bound is 



f(Pn 



2a 



2a + r 



(28) 



For the infeasible case instead, i.e. for < a < follows 
that minimal bound is attained at the boundary point p = -. 
Thus /(2/r) = az(l — r/2) 1 ~ r / 2 . Summarizing, 



2n 



min f(p) 

R3p>max{l,|} 



a > 



2-r 



2a+r — z ^29) 

a% (1 ~r/2) 1 - r / 2 else 



is the best possible upper bound. 

It remains to investigate the conditions for equality. Lieb's 

inequality is fulfilled with equality if strictly p m ; n > - and 
5,7 are matched Gaussians. In this case follows 

r(ax\ + ^x 2 2 ) + {B,x)+C 



(30) 



for some £? £ C 2 , a, C £ C and Re{a} > 0, thus A S7 (-) is 
a two-dimensional Gaussian. Next, to have equality in (16) 
p > 1 and equation (17), which is in this case 

= e Re{-f(oa:?+£a£)+(B,s)+C} 

= Aae -7 <pT7 |:j;|2 = A|C(x)| 7 <F=t 



L S7 (a;)| — e 



(31) 



have to be fulfilled for almost every x. Thus, it follows that 
Re{5} = (0,0), Aae- Re -f c > = 1, Re{a} = 1 and - most 
important - again p = 2^ + 1. But, this is obviously also the 
minimum if a > hence in this and only this case equality 
is achieved. ■ 

It is remarkable that the sharp "if and only if" conclusion for 
Gaussians holds now for a > 2^r. Lieb's inequality alone 
needs a > but in conjuction with Holder's inequality this 
is relaxed. The results are illustrated in Fig.l. Furthermore 
note that for r = 2 every p m [ n is feasible. 
This result is important for so called bosonic Gaussian quan- 
tum channels [9], [12], i.e. C(-) is a two-dimensional Gaus- 
sian. In other words, according to (6) and C(-) as in Corollary 
5, the solution of the Gaussian fidelity problem [13], [3] is 



max Tr{II ff ^(n 7 )} = max X max (A(X)) 



TrJf = l,A'>0 

a 



(32) 



1 



with Gaussian g and 7 as already found in [13] using a 
different approach. But now this states the strong proposition 
that maximum fidelity is achieved only by coherent states. 




Fig. 1. Norm bounds for Gaussian weights: Both functions in (21) separately 
and the combined version are shown for r = 1 and r = 1.9. 



In radar and sonar applications and also for wireless commu- 
nications the following upper bound is important. It is related 
to the case where scattering occurs with constant power in 



some region of phase space (in this context also called time- 
frequency plane). For example in wireless communications 
typically only the maximal dispersions in time and frequency 
(maximum delay spread and maximum Doppler spread) are 
assumed and/or known for some pulse shape optimization. 
Those situations are covered by the following result: 

Corollary 6 Let U c M 2 a Borel set, \U\ < 00 and C(x) = 
mrXu( x ) its C\— normalized characteristic function. Then for 
each p > max{l, |} holds 



rp\U\ 



(33) 



It is not possible to achieve equality. The sharpest bound is 

tS 1 \U\ < 2e/r* 



where r* = max{r, 2}. 



r*|C7| 



r/r* 



else 



(34) 



Proof: The proof is straightforward by observing that 



\C\\ a =\\-^-Xu\\s = \U\- 



According to Thm.4 follows 

f(p)^(-) HCI 

\rp 



rp\U\ 



= e p 



(35) 



(36) 



Equality is not possible because Thm.4 requires C to be 
Gaussian for equality. The optimal version is obtained by 
minimizing the function f(p) under the constraint p > 
max{l, 2/r}. The first derivative /' of / at point p is 



p z rp\U 



)-l) 



(37) 



Thus f (proia) = gives the only point p min = The 
function f(p) i 
h(p)= In f(p) 



function f(p) is log-convex on p £ (0, 2 *\u\ 1 ^ ^' ^ nal ^ s 



-(In ^j^-)/p is convex on /, because 

1 ,. rp\U\ 



ti(p) = -zQn- 



-1) 



^) = i(-21n^ + 3) 
p 6 2 



(38) 



shows, that h"(p) > for all p € I. Hence f(p) is convex on 
/. Obviously p m i n £ I, hence this point is in the convexity 
interval and therefore must be the minimum of /. Further, 
this value is also feasible if still p m i n > max{l, 2/r} = r* jr 
where r* = max{r, 2}, i.e. 

2e 



\U\< 



(39) 



has to be fulfilled. Then the desired result is /(p m i n ) = e™ - ^ . 
If Pmin < r*/r, i.e. is infeasible, the minimum is attained at 
the boundary, i.e. at p — r* jr. Thus 



fir* It) 



r*\U\ 



r fr* 



(40) 



The results are shown in Fig. 2 for r = 1,2,3. For the 
interesting case r = 2 the result further simplifies to 



|A 97 | 



C\\i < 



(41) 



^ \U\<e 
else 

Example: When using the WSSUS model [14] for doubly 




Fig. 2. Norm bounds for j^Xu weights: Both functions in (34) separately 
and the combined version are shown for r = 1, 2, 3. 

dispersive mobile communication channels one typically as- 
sumes time-frequency scattering with shape 



U = {(xi,x 2 ) |0 < xi < 73, \x 2 \ < B d } 



(42) 



with 2B c iJd < 1 < e, where B c i denotes maximum Doppler 
bandwidth Bj and is maximum delay spread. Then (34) 
predicts, that the best (mean) correlation response (r = 2) in 
using filter g at the receiver and 7 at the transmitter is bounded 
above by 



2 C||i <e 



(43) 



III. Conclusions 

In this contribution new bounds on weighted norms of ambi- 
guity functions and Wigner distributions are presented which 
only depend on the shape of the weight function. Further the 
important equality case is discussed which is attained only by 
Gaussian weights and wave functions. The results are impor- 
tant in the field of waveform optimization for non-stationary 



environments as needed for example in WSSUS channels. 
This channel model is frequently used in radar and sonar 
applications and - of course - in wireless communications. 
Furthermore these norms are needed in quantum information 
processing for bosonic quantum channels because they provide 
insights on achievable fidelities in those quantum channels. In 
the special case of the Gaussian quantum channel they provide 
also the optimum input states, i.e. only coherent states achieve 
this optimal fidelity as frequently conjectured. Hence, in the 
mentioned fields the results establish limits on achievable 
performance. 
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